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Abstract 

We discuss lattice methods to obtain the derivatives of a lattice meson mass with respect to the 
bare sea and valence quark masses. Applications are made to quenched and dynamical fermion con- 
figurations. We find evidence for significant differences between quenched and dynamical fermion 
configurations. We discuss how to relate dependence on the bare lattice parameters to more phe- 
nomenologically useful quantities. 



1 Introduction 

In lattice studies of QCD, the action depends on several bare parameters such as the inverse coupling 
(3 and those controlling the quark masses. Here we distinguish the sea quarks which contribute to the 
vacuum and valence quarks which propagate in the vacuum but do not contribute to it. Thus there 
will be two possible parameters describing the quarks masses: the sea and valence hopping parameters 
(k s and k v ). In lattice studies, unlike experiment, it is possible to vary each of these mass parameters 
independently. 

It is of interest to establish the dependence of quantities of physical interest, such as hadron masses, 
on these quark mass parameters. For instance, the valence-quark mass dependence of the meson mass 
controls the J parameter which is related jy to the slope of My versus Mp (where Mp and My arc 
the pseudoscalar and vector meson masses respectively). This slope is found in lattice studies to be 
significantly smaller than the experimental value. It is a challenge for dynamical fermion studies on a 
lattice to narrow this discrepancy as the sea quark mass is reduced. Another area of current interest 
is the magnitude of sea quark effects on hadron masses. The dominant effect of sea quarks is just to 
renormalise the coupling so it is valuable to have techniques to explore in fine detail the sea quark 
effects so that physically significant effects can be explored in dynamical fermion studies. 

One direct way to achieve this is to study the theory at many different combinations of parameters. 
This is the conventional way to study the valence quark mass dependence and is reasonably efficient 
since the lattice configurations themselves do not depend on n v . For the sea quark mass, however, 
this is a computationally challenging endeavour since different gauge configurations must be constructed 
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for each k s value and then the finite differences of hadron masses between these different ensembles of 
configurations will be small and quite noisy. 

One way to obtain estimates of derivatives by working with a lattice ensemble at one set of parameters 
is described in refQ. Here we specialise to explore a method to obtain the derivative of a hadron mass 
with respect to a parameter such as k s . The method is essentially to take formally the derivative of a 
lattice identity. This method, often called a 'sum rule', has been used before to obtain derivatives with 
respect to j3 ||. Here we use a similar approach to extract derivatives with respect to k s and k v - see 
also J§. 

The derivative with respect to k v involves a three point function of fermion fields and so cannot be 
obtained from propagators from one source only. Here we choose to use stochastic propagators |(| with 
maximal variance reduction Q which allow the appropriate propagator combination to be evaluated. 

For the derivative with respect to k s , a disconnected three point function is needed. In this case we 
use Z2 noise methods J?], || to evaluate the appropriate combination of propagators. We apply this to 
quenched and dynamical fermion gauge configurations and see a significant difference. We discuss the 
impact of these results on the sea-quark dependence of meson masses. 

This study is exploratory and we discuss the computational effort needed to extract these derivatives 
with respect to bare quark masses. We also compare our results with those obtained by taking finite 
differences. 



2 Quark mass dependences 

The mesonic masses in lattice studies are determined by measuring two-point correlations of appropriate 
operators at large time separation t. We then wish to take the formal derivative with respect to a 
parameter representing the quark mass. This will give the required sum rules for the derivative of the 
lattice hadron mass with respect to the quark mass parameter. 

Consider an action density 

N f 

S=Sf+ l3Sg = J2 + PSg (I) 

1 

where, for the Wilson-Dirac discretisation of fermions, 

M = m + D (2) 

where the quark mass parameter m = with k the conventional hopping parameter, so in terms of 
the bare quark mass mj, in the naive continuum limit m = 8 + 2arrib- The term D contains the Wilson 
nearest neighbour gauge link terms as well as the SW-clover terms with coefficient Csw- The hadronic 
correlation is then given by 

C{t) = I J H(0)H^t)e s (3) 

Here for mesons H will be of the form tpTip and in this work we will concentrate on the case of flavour non- 
singlet mesons so that the hairpin diagrams will not be needed. Then the fermionic degrees of freedom 
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are integrated out giving a factor of the inverse of the fermion matrix (G = M 1 ) for each pairing: 

C(t) = (0\(G(0,t)TG(t,0)T)\0) (4) 

At large t, this correlation will be dominated by the ground state meson with the quantum numbers 
created by H: 

C(t) = c 2 e- Mot + ... (5) 

This sketch of the formalism allows us to explore taking the derivative with respect to the quark mass 
parameter m (actually the inverse hopping parameter) on each side of the above expressions for C(t). 
This derivative is to be taken at fixed (3. Then, since formally the only m-dependence is in the exponent, 
the derivative brings down a factor of Nftptp. 



dm 



= ± J H{Q)H\t)N f ^e s - C(t)l J N f ^e s (6) 

where the second term comes from the m-dependence implicit in Z. On integrating out the six fermions, 
this will give two diagrams, connected and disconnected (actually only the connected part of the discon- 
nected diagram will contribute as discussed below). Thus, summing explicitly over the insertion at t\, 
we have 

^ = E (~ C 3 (*i . t) + N f D 3 ( tl ,t)) (7) 
ti 

o 

t, tOt! t 

/ X_ \ / \ 

r r r r 

y \ / 



(a) Connected Mesonic Correlator (b) Disconnected Mesonic Correlator 

Figure 1: 

The diagrams corresponding to C3 and D 3 . 

Where for the connected diagram there will be terms from the insertion on either quark line: 

C 3 (t u t) = (0| (G(0,ti)G(ti,t)TG(t,0)T) + (G(0,t)rG(t,ti)G(ti,0)r) |0) (8) 
while the disconnected diagram is, for each flavour of quark in the loop, 

D 3 (h,t) = (0| (G(0,t)TG(t,0)T) (G(h,h)) |0) - (0| (G(0,t)TG(t,0)T) |0)(0| (G(ti,*i)) |0) (9) 
where curved brackets imply a trace over colour, spin and space coordinates. 
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In the quenched approximation, only the connected diagram contributes. This can be seen another 
way since M.G = 1 implies M.dGjdm + dM.jdmG — and since dAi/dm = 1 by definition, then 
dG/dm = —Mr x G = —GG. Thus either fermion propagator in the mesonic correlator can be 'opened' 
by an insertion. 

For dynamical fermions, both types of diagram contribute but one can see that the connected diagram 
corresponds to varying the valence quark mass while the disconnected diagram corresponds to varying 
the sea quark mass. 

The two point hadronic correlation can be expressed in terms of a sum over intermediate states of 
masses Mi 

C(t) ="£c*e- Mit (10) 

i 

The leading term in the derivative at large t can be then be evaluated 

dC{ t) = _£^V e -M t + d£o e -M t + _ (n) 
dm dm dm 

and it thus behaves as te~ Mot where Mo is the ground state meson mass. 

We now extract the contribution from the right hand side which has this same behaviour. For both 
C3 and D3, the insertion is summed over all space and time. Then the leading term arises when the 
lightest allowed meson propagates and when < t\ < t. This will produce terms which are linear in 
t which arise from the t possible insertions (at t{) between the creation and destruction of the meson. 
Then evaluating this ground state meson contribution, for the connected diagram, gives 



E 

ti 



C 3 (h,t) = c 2 t(x^+X^ e~ Mot = t (iC'+I 12 ') C(t) (12) 



where the suffix refers to the insertion on quark propagator 1 or 2 and X is the matrix element of the 
iptp insertion between ground state hadrons. A similar expression applies for the disconnected case. 



Equating the coefficients of the terms behaving as te M °* on each side of the identity, then gives the 
exact result that 

,1 1 r 

iW + lM (13) 



dM 



dm v 

where the matrix element sum can be obtained by extracting the ground state contribution to C3/C. In 
principle this can be obtained by taking the insertion in C3 such that < t\ < t and both t\ and t — t\ 
are large so that the ground state contributes. So we can write 

^= Bm (14) 
dm v ti,(t-*i)— >oo C(r) 



This sum rule relates the derivative to an expression that can be evaluated from lattice configurations 
at only one set of parameters. It is an exact identity. If there is a dependence of the lattice meson mass 
Mo on the finite spatial size L of the lattice, the derivative should be taken at fixed number of lattice 
spacings, not at fixed physical size. These considerations are very similar to those used in the lattice sum 
rules derived by taking formal derivatives with respect to [3 S. 
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For the disconnected diagram, the equivalent expression is 



d -^ = -N f Inn ( 15 ) 

dm s ti,(t-ti)-»oo C(t) 



In order to evaluate these expressions on a lattice, it is useful to consider efficient ways in which 
excited state contributions can be eliminated, since the formal limits of large t will have big noise to 
signal. Here we consider the connected correlation C3 and use a complete set of hadron states of mass 
Mi in the intermediate intervals of time extent t\ and t 2 = t — t\. 

In practice, we will be using more than one operator to create and destroy the hadronic state. This 
allows an optimal combination of these operators to be formed that minimises the excited state contri- 
bution. Then the two-body correlation between operators a at t = and b at t will be given by 

C^ b \t) = Y,c < f ) e- Mit cf ) (16) 

i 

Cf b) (t u t 2 ) = J2^ ) e- Mitl x ij e- M ^cf ) (17) 

where Xoo is the required quantity (X^ + X^) - the matrix element appropriate to the ground state 
meson of mass M . We might expect Xn to be similar in sign and magnitude to x 00 if the quark 
mass dependence of the excited state is comparable to that of the ground state and thus excited state 
contributions would cancel in the ratio C3/C. This is incorrect, since the off-diagonal terms (xoi) will 
dominate the excited state contributions to C3 since the excited state only propagates for the shorter 
interval t\ (or t 2 ). One way to extract xoo is to make a fit to the three point data with both t\ > t m ; n and 
t2 > £min, keeping the coefficients and masses (cf^ and Mi) fixed from the fit to the two-point function 
data with t > i m i n . This can be compared with the more direct approach of looking for a plateau in 
C3(ti,t 2 )/C(t) as t\ and t 2 are increased (with t = t\ + t 2 ). 

When two (or more) different types of hadronic creation operators are used, a variational method is 
an effective way to determine the ground state contribution to C and hence to extract the ground state 
contribution to C3 . Alternatively, if a two state fit to the two-body correlation between two operators at 
each end is made, then from the coefficients, it follows that a combination of operators c^ 1 H\ — H 2 
will remove the contribution of the excited state in the approximation that only two states contribute to 
the correlations. Then we can use this combination to evaluate the ground state component of C3/C, 
using 

(2) (2)^(11) _ „ (2) (1W12) (1) (1)^(22) 

c (2) c (2) C (iD - 2cl 2) 4 1) C7(i2) + c^c^CW 1 ! 

A similar analysis holds equivalently for the extraction of the ground state contribution efoo to D 3 . 

When t i=a T/2, the contributions from propagation around the time boundary of the lattice may be 
significant. For C3 there will be no such 'round the back' term because the insertion is made explicitly, 
while for D3 the connected matrix element involved will cancel for the round the back term. In contrast 
the two-body correlator C will be a sum of two terms. Illustrating this for the ground state component 
for one type of operator, we have: 

C - cge- Mot + cge-^-O (19) 
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G 3 = c 2 e- M ^x 00 (20) 

Hence 

C3/C = SM 1 + e J 0(r _ 2t) (21) 

This formalism can be used to correct for the different i-dependences when looking for a plateau in G3 jC 
and in D3/C as t increases. 

We now discuss efficient methods to evaluate these correlators on a lattice. 



3 Valence quark mass dependence 

As a first application, we consider the dependence of the hadron mass on the valence quarks. For this a 
three point function needs to be evaluated - see fig. fa. Thus conventional quark propagators from one 
source are inadequate for this task. One feasible way forward is to use a stochastic inversion method 
which allows the evaluation of quark propagators from any site to any other site. Although the stochastic 
method is not more efficient than the conventional inversion from one source for mesons made of light 
quarks it does allow the flexibility to evaluate three point correlations readily. For this reason it 
allows an exploratory study of this area. 

Stochastic propagators |5|, ^| are one technique to invert the fermionic matrix for the light quarks. 
They can be used in place of light quark propagators calculated with the usual deterministic algorithm. 
The stochastic inversion is based on the relation: 

Gh = MTj 1 = \\t>4> {M jk <t>k)*^i exp (-0* (M^M)ij(pj) (22) 

where, in our case, M. is the improved Wilson-Dirac fermionic operator and the indices i,j, k represent 
simultaneously the space-time coordinates, the spinor and colour indices. For every gauge configuration, 
an ensemble of independent fields <f>i (we use 24 following jq|) is generated with gaussian probability: 

Pffl = | exp (-tiiM^M)^) (23) 
All light propagators are computed as averages over the pseudo-fermionic samples: 

Ga = { or (24) 
{ l5 {<f>*{M<t>)i)^ 

where the two expressions are related by Gy = 756^75. Moreover, the maximal variance reduction 
method is applied in order to minimise the statistical noise Q . The maximal variance reduction method 
involves dividing the lattice into two boxes (0 < t < T/2 and T/2 < t < T) and solving the equation of 
motion numerically within each box, keeping the pseudo-fermion field 4> on the boundary fixed. According 
to the maximal reduction method, the fields which enter the correlation functions must be cither the 
original fields <j) or solutions of the equation of motion in disconnected regions. The stochastic propagator 
is therefore defined from each point in one box to every point in the other box or on the boundary. For 
this reason, when computing the three-point correlation function, 

(0\H(h, x)O(t , y)H^(t 2 , z)\0) (25) 

x,y,z 
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the operator O (which is {flip) is forced to be on the boundary (to — or T/2) and the other two 
operators must be in different boxes, while the spatial coordinates are not constrained. If j is a point 
of the boundary, not all the terms in (Mcj))j lie on the boundary because the operator M. involves first 
neighbours in all directions. Hence, whenever a propagator Gij is needed with one of the points on the 
boundary, we use whichever of the two expressions in Eq. ^4] has M.c\> computed away from the boundary. 
This implies that we are restricted to t > 2. 

The numerical analysis used 24 stochastic samples on each of 20 quenched gauge configurations, 
generated [|| on a 12 3 x 24 lattice at (3 = 5.7, corresponding to a -1 = 0.91 GeV. With improved clover 
coefficient Csw — 1-57, we use two values of n: n\ = 0.14077 and K2 — 0.13843. The lighter value K\ 
corresponds to a bare mass of the light quark around the strange mass. The chiral limit corresponds to 
k c = 0.14351 0. Error estimates come from bootstrap over the gauge configurations. We also made an 
exploratory study of some dynamical fermion configurations, as will be discussed later. 

In smearing the hadronic interpolating operators, spatial fuzzed links are used. Following the pre- 
scription in || [To), to which the interested reader should refer for details, the fuzzed links are defined 
iteratively as: 

[/new = V (fU old + Ubend,ij (26) 

where V is a projector over SU (3), and f/bend.i are the staples attached to the link in the spatial directions. 
Five iterations of fuzzing with / = 2.5 are used and then the fuzzed links are combined to straight paths 
of length three. The fuzzed fermionic fields are defined following |Io[ . 

We employed two types of hadronic operator for the correlations - local and fuzzed - yielding a 2 x 2 
matrix. From this we use a variational approach to extract the linear combination of operators which 
maximises the ground state contribution - as described above. Since we are able to get good two state 
fits to the two-body correlations for the pseudoscalar meson for t > 3, this variational linear combination 
was determined using i-values 3 and 4. In order to maximise the ground state contribution relative to 
excited states, we evaluated the three point diagram Ca(ti, t-z) using values of t\ and ti near to t/2 where 
t = t\ + ti. The ground state improved ratio of C§(t\, t2/C(t) is plotted in fig. ^. The extraction of the 
ground state should be good if ti,t 2 > 3. For odd values of t there are higher statistics (from the 3,4 
and 4,3 partitions of t=7 for instance). Thus we expect t — 7 to be the best determined value and this 
is given in table |l|. Consistency at higher i-values confirms that the ground state extraction is correct. 

Table 1: Connected loop correlations. 

K v Mp dMp / dm v My dMyjdm v n gaugc 



0.14077 0.529(2) 1.97(27) 0.815(5) 1.9(9) 20 

0.13843 0.736(2) 1.56(22) 0.938(3) 0.3(4) 20_ 

0.1395 0.558(8) 1.3(3) 0.786(9) 0.2(1.1) 5 



For the pseudoscalar meson, we expect that Mp is approximately linear in m v . Thus dM p /dm v 
should decrease like 1/Mp which is indeed consistent with the results shown in table [l]. From the 
high statistics spectroscopy at these two hopping parameters [pL one can evaluate the finite difference 
obtaining dMp/dm v = 2.18(4) which agrees very well with the values determined from the sum rules of 
2MpdM P /dm v = 2.08(29) and 2.30(32) at k v = 0.14077 and 0.13843 respectively. 
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Figure 2: The connected correlation C3/C versus time t = t\ + t% in lattice units. The data are for the 
variational combination that reduces the excited state contribution and are from quenched lattices with 
k = 0.14077 with \ti —ti\ < 2. We expect the ground state contribution to be dominant when t\ > 2 and 
t% > 2, that is for t > 6. 



For the vector meson, the expectation is that My is approximately linear in m v . The finite differ- 
ence [|| gives dMv/dm v = 1.02(7). The sum rule determination with our current statistics is too noisy 
at t = 7 to give an accurate value. For the heavier quark mass (k.„ = 0.13843) a good two-state fit 
to the two point correlation data can be made for t > 2. This allows us to use t = 4 and 5 for the 
C3/C ratio and we obtain 0.90(18) and 0.93(26) respectively, in excellent agreement with the expected 
value. At the lighter quark mass, we need t > 3 for a two-state fit so the poor result remains. This is 
a disappointment, since from the values of dMy /dm v and dMp/dm v , one can evaluate the J parameter 
(which is the physical quantity, defined in the continuum limit as MydMy j 'dMp at My/Mp = 1.8) at a 
quark mass corresponding to m v . Thus J can be determined at the lightest quark mass directly, rather 
than as a difference between two quark masses. This J parameter is a useful indicator [0 of the distance 
between quenched QCD (with J « 0.37) and experiment (with J = 0.48(2)). Hence a quick and accurate 
method to determine J would be useful to calibrate dynamical fermion studies. 

We also evaluated the same quantities for dynamical fermion configurations |ll[] at (3 — 5.2 with 
two flavours of sea quarks at k s = 0.1395 on a 12 3 24 lattice using a SW-clover improved action with 
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Csw — 1-76. The correlation was evaluated with k v = k s . In this case the higher statistics determination 
of the masses allows the derivative at fixed k s to be evaluated, giving dMp/dm v = 2.2(4) from the 
finite difference between k v of 0.1395 and 0.1390. Our analysis is from only 5 gauge configurations and 
so the error may be underestimated because of the small sample size. For the pseudoscalar, we find 
acceptable two state fits for t > 2 (the hadronic operators are local and fuzzed with straight paths of 2 
links) and so we may use the variational method from t of 2 to 3 to determine the ground state couplings. 
This yields a value at t=4 and 5 of C3/C of 1-4(4) and 1.5(3), respectively. This is consistent with the 
value in table |l] and with the finite difference value within errors. The vector meson case is too noisy to 
be of any use. The main conclusion is that the ratio of correlations C3/C is very similar in the dynamical 
configurations to the quenched case. This is not really surprising since the sea quark masses used in the 
dynamical quark study are fairly large - larger than the strange quark mass. 



4 Sea quark mass dependence 

The disconnected diagram (see fig. lb) involves measuring two gauge invariant contributions: the two- 
point hadronic correlator C and the loop contribution corresponding to TraccM -1 where the trace is a 
sum over colour, spin and spatial coordinates at a given time to- This needs the propagator from each 
site on a time slice to a sink corresponding to the same site. There is an efficient way to evaluate this 
making use of Z2 stochastic sources |?], |§| . Here we propose a variant of this method which is appropriate 
for our current study. This method also gives the two point correlator C{t\ 1 t2) for pseudoscalar mesons 
and vector mesons from any time t\ to any other time ti. Then combined with the loop contribution 
at to, we have the ingredients needed to evaluate the required connected part -D3 of the disconnected 
correlation. 

Details of the Z2 method used are given in the Appendix. In this exploratory study on 12 3 24 lattices, 
we use local operators to create the pseudoscalar and vector mesons. We have used rather generous 
values of the number of Z2 samples per time slice (namely between 16 and 32 for each of the two related 
types of source used) . This amounts to 768 or more inversions (equivalent to 64 conventional propagator 
inversions from 12 colour spin sources) per gauge configuration. Because of the decreased number of 
iterations of the inversion algorithm in our case, the time used is equivalent to about 30 conventional 
propagator determinations per gauge configuration. This is a substantial computational challenge, but 
it does provide a significant resource: the loop contributions at each t and the pseudoscalar and vector 
correlators from any t\ to any ti. Because of our choice of number of Z2 samples, we have negligible 
errors coming from the Z2 noise for the value of Traced -1 from each time-slice and for the pseudoscalar 
correlator from t\ to t^- For the vector meson correlator, the error from the Z2 method is in some cases 
comparable to the intrinsic variation and we correct for this in derived quantities by increasing our errors 
appropriately where necessary. Indeed, in retrospect, it would have been more efficient for the present 
study to use less Z2 samples and to explore more gauge configurations. Our approach, however, was that 
so much computational effort has gone into the production of the dynamical fermion configurations that 
the large number of inversions used in measurement are in effect a relatively small extra overhead. 

We evaluated these quantities for dynamical fermion configurations JllJ at ft = 5.2 with two flavours 
of sea quarks at n s = 0.1390, 0.1395 and 0.1398 on a 12 3 24 lattice using a SW-clover improved action 
with Csw — 1.76. In our evaluations we restrict ourselves to the case where the propagating quarks have 
the sea-quark mass, i.e. k v = K s . The number of gauge configurations used and number of Z2 samples 
nz are given in table ||. We also quote, for completeness, the pseudoscalar and vector meson masses and 
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the Rq values obtained from higher statistics by conventional methods |j| [L2| [lT| . 

It is possible to measure the disconnected diagram in quenched gauge configurations as well as in 
dynamical fermion configurations. We use the same gauge configurations as discussed in the previous 
section. For the quenched case, we include a factor of Nf = 2 explicitly to facilitate comparison with the 
dynamical fermion configurations that have Nf = 2. 

As discussed previously, the disconnected 3-point correlation D^{tx,t2) can be fitted to obtain the 
matrix element doo that gives us dM/dm s . Because we only have data on the correlations from local 
hadronic operators in this study, we choose to make use of the results of conventional studies of the 
2-point correlators from both local and non-local (smeared or fuzzed) operators from larger samples of 
configurations |llj to determine the couplings Cj of the ground state and excited state mesons to our 
operators. We find that adequate two-state fits can be made to these 2-point correlations for t > 2. Then 
keeping the masses and coefficients Cj fixed, we can fit all the 3-point data with t\ > 2 and £2 > 2. Some 
typical fits are shown in fig. |^. The fit results are shown in table |J the upper two lines are from quenched 
configurations while the lower three lines are with dynamical fermions. 



Table 2: Disconnected loop correlations. 



K 


M P 


dMp 1 dm s 


My 


dMy j dm s 


i?0 


^gaugc 


n z 


0.14077 


0.529(2) 


1.18(26) 


0.815(5) 


1.6(5) 


2.92(1) 


20 


16 x 2 


0.13843 


0.736(2) 


0.86(15) 


0.938(3) 


1.0(2) 


2.92(1) 


20 


16 x 2 


0.1398 


0.476(14) 


3.0(5) 


0.706(16) 


2.7(6) 


3.65(4) 


20 


32 x 2 


0.1395 


0.558(8) 


3.1(6) 


0.786(9) 


3.0(7) 


3.44(6) 


20 


32 x 2 


0.1390 


0.707(5) 


1.9(4) 


0.901(10) 


1.8(3) 


3.05(7) 


24 


20 x 2 



The sign of the effect implies that the loop (T = TraccM 1 ) is anti-correlated with the pion two-point 
correlation C which straddles it in time on a lattice. This anticorrelation is large with, for example, 

(SC ST)/(((SC) 2 } {{ST) 2 )) 1 ' 2 « -0.5 (27) 

at t = 6 for both the dynamical fermion and quenched cases. 

This anti-correlation is seen to be very similar for pseudoscalar and vector mesons. One qualitative 
argument for the sign of the correlation is that an upward fluctuation of C corresponds to configurations 
in which quarks propagate easily over large distances whereas an upward fluctuation of the loop (T) 
comes from configurations in which quarks do not propagate easily - and so have a bigger amplitude 
at the origin. In terms of our identities which relate this disconnected correlation to the derivatives 
dM/dm s , we see that the main effect comes from the dependence of the lattice spacing a on m s at fixed 
(3. It is well known that a decreases as the sea-quark mass decreases: indeed this is why the (3 value 
used in dynamical simulations is smaller than that used in quenched. The UKQCD study jll| of the 
dynamical fermion configurations we are using finds d \oga/dm s « —4 - as shown in fig. |J. Furthermore 
the slope appears larger at smaller sea quark mass - in line with what we find in table |[ 

We also measure the same disconnected correlation in quenched configurations. The results are 
qualitatively similar to those from dynamical fermion configurations. This implies that one can explore 
the sea-quark dependence of meson masses using quenched configurations. This appears a striking advance 
- one can get at essential information concerning sea-quarks without the heavy computational overhead 
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Figure 3: The disconnected correlation —NfD^/C with Nf = 2 versus time t = t\ + £2 in lattice units. 
The upper data are from dynamical fcrmions with K soa = 0.139, while the lower data from quenched 
lattices with n — 0.13843. The curves show the two-state fits to these data with \ti — £2! =0 or 1 as 
described in the text. The additional points (crosses and octagons) have \t\ — t2 \ = 2 or 3 and are fitted 
by the dotted curve. 

of dynamical fermion simulations. However, it is widely appreciated that most lattice obscrvables are 
insensitive to the presence of sea quarks if the lattice spacing and MpjMy ratio are lined up. Hence, 
once one has expressed the quantity of interest as a vacuum expectation value, it may be evaluated using 
quenched configurations. We now explore this in a little more detail. 

For the heavier quenched (k = 0.13843) and dynamical (n — 0.1390) cases, the lattice spacings (taken 
from Rq 3) and the MpjMy ratio (at 0.78) are very similar. Thus we may directly compare the 
dM/dm s values obtained. From table |2j, we see that dM/dm s has significantly smaller values (by two 
standard deviations) in quenched than in dynamical fermion configurations. This conclusion is reinforced 
by the presentation of the fits to these data shown in fig. |^. These data suggest that this observable 
is indeed capable of distinguishing between configurations with different sea quark structure. One note 
of caution is that since the lattice spacing is rather coarse, the different finite lattice spacing effects in 
quenched and dynamical simulations may be partly responsible for this observed difference. 
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Figure 4: The lattice evaluation of Rq for dynamical fcrmion configurations with sea quarks of hopping 
parameter k from rcfjll|]. We define m s = 1/re. 

This ability to distinguish quantitatively between quenched and dynamical gauge configurations is 
important - in most cases previously studied, no such discrimination was detectable. That the observable 
currently under study allows this discrimination is not entirely unexpected since the dynamical fcrmion 
configurations are weighted by det(A4) which is closely related to TraceVW -1 which is a component of 
the disconnected correlator. 

For dynamical fermions it is possible to evaluate by conventional methods the hadronic mass differences 
as the sea quark mass is varied and so obtain an estimate of the derivative which can be compared with 
our results. For the pseudoscalar meson, finite difference determinations |fjl| at fixed valence quark mass 
of k v = 0.1390 give dM P /dm s = 3.1(3) and of n v = 0.1395 give dM P /dm s = 4.3(4). Both of these 
finite difference estimates are somewhat larger than the derivatives determined above. The situation is 
the same for the vector meson mass derivatives where the finite difference determinations give 3.4(5) and 
4.0(1.2) at a fixed valence mass of 0.139 and 0.1395, respectively. 

The two different approaches to determining these sea-quark mass derivatives used different gauge 
ensembles (propagators from the origin from about 100 gauge configurations for ref ill) , compared to 
propagators from all sites on about 20 gauge configurations here) and the differences are only at the two 
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standard deviation level. At present the quoted statistical errors from the derivative method we use are 
comparable to those from finite differences of masses. Using the full set of gauge configurations available, 
our derivative method would give the more accurate determination of the sea quark dependence of the 
meson masses. 



5 Discussion 

There are several issues of interest in determining the dependence of hadron masses on the quark masses. 
Here we are not concerned with the problem of defining precisely the quark masses. Rather we discuss 
the dependence of the meson masses as the sea quark mass is reduced to look for explicit signs of different 
physics as the quark loops become more important in the vacuum. One of the complicating features in 
the lattice approach is that changing the sea-quark mass parameter has several consequences - among 
them that the lattice spacing is changed. 

As an illustration, since the lattice spacing a depends on the sea quark parameter m s , let us consider 
the dimensionless ratio My/Mp. Then 

dm s n( M P ' My dm s M P dm s 

can be evaluated. Since we find dMp / dm s ss dMy j 'dm s , this gives a negative result which implies that 
the My/Mp ratio increases as the sea quark mass is decreased. This change in sea-quark mass parameter 
is at a constant n v , however, which is not necessarily what is required. 

To clarify this discussion, it must be remembered that the bare parameters (/3, k v ,k s ) which occur in 
the lattice formalism are not simply related to the more physical parameters a and the sea and valence 
quark masses which we denote here as q v and q s . One example of this intricate relationship is that as 
n s is increased (i.e. towards n c so that the sea quarks are lighter) then a becomes smaller (this can be 
seen from the observation that (3 needs to be reduced for dynamical fcrmions to keep a approximately 
the same) . Furthermore, this change of n s is also likely to result in a different value of n c (here defined as 
the value which gives a massless pion on varying k v at that sea quark mass) and hence the relationship 
of k v with q v will be modified too. 

Thus one needs to set up a prescription to determine appropriate values of the lattice parameters. 
One proposal is to identify physical quantities which should not depend on all of the physical parameters. 
Thus we can choose to use ro (defined via the static potential at moderate separations) to determine the 
lattice spacing a, assuming it to be independent of the quark masses. This would not have been true if 
the string tension were to have been used to set the scale since the string breaking at large separation 
will be strongly affected by the sea quark mass. 

For the quark mass dependence, we arc considering a world where the valence quark mass q v can be 
varied independently of the sea quark mass q s . This is not so far from experiment if one regards the u : d 
quarks as sea quarks and the strange quark as a valence quark whose contribution to the sea is relatively 
small. 

To isolate the quark mass dependence, we choose to make use of a very conspicuous experimental 
fact: the vector mesons are 'magically mixed' with the 4> meson being almost pure ss while the p and lu 
are almost degenerate and composed of u, d quarks. Furthermore the <fr has much reduced decay matrix 
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elements to final states containing only u, d quarks. This is the OZI rule: disconnected quark diagrams 
are suppressed. All of this phenomenology suggests that the vector meson nonet is well described by 
the naive quark model: it does not contain significant sea quark contributions to the masses. Thus we 
choose to define the sea quark mass q s such that the vector meson masses are independent of it. For 
other mesons, especially the pseudoscalar mesons, we do expect some dependence of the masses explicitly 
on the sea-quark mass q s and we shall try to estimate it. 

One way to proceed is to remove the explicit a-dependence of the lattice masses by forming the product 
with Rq. Then RqMp will be equal to the continuum product romp up to lattice artefact corrections 
which are of order a for the Wilson fermion discretisation but the clover-improvement scheme we use 
should reduce these lattice artefact corrections to being dominantly of order a 2 . For ease of notation we 
define P s — d(R Mp)/dm s etc. Here we assume, as discussed above, that i?o is independent of m v and 
that it does depend on m s through the dependence of a on m s . This sea quark mass dependence of Rq 
can be extracted by explicitly evaluating Rq at a range of m s values p"i"| ] , as illustrated in fig. ^, giving 

' " /f> ° = -4.7(1.8) ; -4.0(2.0) (29) 



Rq dm s 



where differences are taken from n s of 0.1390 to 0.1395 and then 0.1395 to 0.1398 respectively. These 
values can then be used to obtain 

1 1 d{R M P ) _ 1 dR 1 dM P 



RqMp RqMp dm s Rq dm s Mp dm s 

where a substantial cancellation occurs between the latter two terms. Thus we find that the resulting 
errors are sufficiently large that even the sign of P s is not well determined. However, the sign of P s does 
not necessarily have any direct physical meaning as we now discuss. 

Assuming one had accurate values, we now discuss how to interpret them. The situation is illustrated 
on a plot of P = RqMp against V = RqMv in fig. |[ 

As the valence quark mass parameter m v is varied, a curve is traced out. What is of interest, however, 
is the difference between such curves as the sea quark mass parameter m s is varied. Assuming, as discussed 
above, that the vector mass (V = RqMv) is independent of q s , then yields the required dependence of 
the pseudoscalar mass on m s at fixed V: 



d{R M P ) 



dm s 



Ps - JT-Pv (31) 



We claim that this quantity will give the physically relevant part of the sea-quark dependence of meson 
masses: d(RoMp)/dq s — r d(mp)/dq s ). Indeed a presentation in this spirit was already shown in ref. |ll) . 
There it was concluded that as the sea-quark mass is reduced, the meson masses move towards closer 
agreement with the experimental data point (t^ s , <fi) with d(R()Mp) /dm s \v > (here rj s is the mass 
expected for a ss pseudoscalar meson). 

Since the precision we obtain in this preliminary study on the derivatives is not superior to that 
which was obtained by directly varying the sea and valence masses ]Tl"| , the conclusions of that work are 
not modified. However, for dynamical fermion studies where only one sea-quark mass is employed, our 
methods will enable the derivatives with respect to the sea-quark mass to be evaluated. 
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V(=RoM v ) 

A 



( P v> V v) 




P(=R<M,) 

Figure 5: An illustration of the bare valence (v) and sea (s) quark mass dependence of pseudoscalar (P) 
and vector (V) meson masses in units of i?o- 

6 Conclusions 

One of the current problems in lattice study of hadron spectra, is to evaluate the physical consequences 
of including sea quark effects in the vacuum. We have presented lattice techniques to evaluate the 
dependence of meson masses on the valence and sea-quark parameters. These techniques allow such 
studies to be made using gauge configurations at a single set of lattice parameters. This is a significant 
advance for dynamical fermion studies which are very computationally intensive. Moreover it implies 
that some estimates of these sea-quark properties can even be made using quenched configurations. 

We have discussed how to extract physically useful information about the sea-quark effects from these 
observables. Our proposal takes account of the changes induced in the lattice spacing and in the valence 
mass definition as the sea quark parameter is changed. 

One rather encouraging feature is that we see evidence for a significant difference for the disconnected 
correlation ratio (our D3/C) between quenched and dynamical quark configurations. It will be of interest 
to explore this difference at finer lattice spacing to establish that it is indeed a continuum effect. 
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A Appendix: Analysis of Z2 noise vector methods. 



A.l Introduction 

We summarise first the salient ideas in the Z2 method (7[ ||, before indicating the special features that 
we have made use of. 

The required time-slice loop term can be expressed in terms of the quark propagator .M -1 on a given 
gauge configuration as 

T(t) = Y / M-l aj (x,t;x,t) (32) 

X 

where we explicitly show the colour index a and Dirac index j here. Since M. is 75-hermitian, then T 
is real on any time-slice of any gauge configuration. To evaluate this expression for all a; on a time slice 
using point sources would require solving the lattice Dirac equation for I? sources of each colour and 
Dirac index. Let us instead explore using distributed sources £ p (:r,i)bfc where p labels the source. 

Then solving the lattice Dirac equation from such a source 

G p aj (x',t') = M-l bk (x',t';x,t) e bk (x,t) (33) 

and combining with an appropriate combination involving the same source, we have 

r^)^^( Il trG;.( I , i )^;.( Il tr^ tt (M;x' 1 t') e bk (x',f) (34) 

X 

Interpreting the sources £ as random with specific properties then makes this quantity, averaged over 
realisations of the random source, to be just that required, namely 

(7*(t)) = Tit) (35) 

This allows the possibility of an unbiased estimate of the required quantity with a moderate number of 
inversions inz)- We require that the random sources (,^Jx,t) with p = 1, . . .nz are such that the only 
non-zero expectation values of bilinears are given by 

(Cai^l)*!)* £bfc(£2,t2)) ^ S pq d ab S jk 6 Xl . X2 d tl . t2 (36) 

This can be implemented by assigning an independent random number to each site x, colour a and Dirac 
index j for each sample p. The optimum distribution of those random numbers can be chosen to minimise 
the variance of the required observable. 

The variance of this estimator is minimised @ by taking Z2 noise (more correctly Z2 x Z2) , namely 
each component (for real and imaginary parts separately) to be randomly ±1/^/2. Then 

a\ = iRcal (.V1, ; '.V-1, : + M^MJ*)^. (37) 

where only the off-diagonal part of contributes and here we include space, time, colour and Dirac 

indices into i. 

The variance can be reduced by using a more selective source, for example ^| with specific Dirac 
components. This involves more inversions, however, if the full signal is to be evaluated. Here we choose, 
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instead, to use a source which is only on a specific time-plane to- Thus in the above formalism £, p {x,t)bk 
is to be taken as zero outside the time-slice to of interest. This reduces the variance by a factor of 
approximately 4 at the expense of 24 (in our case) times as many inversions. This is not cost-effective 
for evaluating T but it does enable us to extract mesonic two-point correlators as we now discuss. 

A. 2 Meson correlators 

It is also possible to use Z2 source methods to determine meson correlators. For illustration, consider 
the correlator between local hadron operators of zero momentum given by the average in the gauge 
configurations: 

C(t) = (0|#(ii)tft(t 2 )| ) (38) 

with t = \ti — *2 1 and where 

H(t) = J2 (x, t)T jk i> ak (x, t) (39) 

X 

creates a meson with quantum number given by the Dirac matrix T, where r = 75 for pseudoscalar 
mesons and V = ji for vector mesons. 

Then, using the 75-hermitian property of the fermion matrix, we need to evaluate (suppressing the 
colour indices) 

C(t) = (0|(T7fe)«Af«(a:i,ti;a:2,t2)(r'ys)fci J M^ 1 *(a;i ) ti;a:2,t2)|0) (40) 

Xl,X2 

This can be evaluated using Z2 methods where G p (t) is the propagator from source £ p on time-slice 
t provided one also has the propagator G pT (t) from source (r7 5 )£ p on the same time-slice. Then the 
average over samples p of this source will give the contribution to C(t) from one time-slice on one gauge 
configuration: 

G{t) - (52G p i (t 1 -,x 2 ,t 2 )Gf*(tv,x 2 ,t 2 )(T l5 ) ij ) (41) 

This method allows us to obtain mesonic correlators from any time slice to any other. For the pseudoscalar 
meson, no additional inversions are needed since 1^75 = 1 in that case. For the vector meson case, we use 
r = 7i with i = 1, 2 or 3 randomly chosen for each sample p. 

In principle, one could obtain mesonic correlations using Z2 methods without additional inversions - 
for example by explicitly evaluating the average over samples p, q of 

Y,Gi(ti;x2,t2)T ij e j *(x2,h) x^Gl^lxut^Tk^ixuh) (42) 

X 2 Xl 

In this case, combinatorial factors make the variance of this estimator comparable to the signal so it is 
an inefficient estimator. Using sources at all t would aggravate this problem considerably. 

As described in the text, we can combine the Z2 estimate of the loop at to with the Z2 estimate 
of the mesonic correlator, provided to is not the source point of the mesonic correlator determination. 
This restriction is of no consequence since we are interested in a loop roughly midway along the mesonic 
correlator. 
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A. 3 Propagators from Z2 sources 



The techniques used to evaluate the propagator from a given Z2 source are just those used in a standard 
inversion from any source. This is achieved by an iterative inversion process (either minimal residual or 
BiCGstab algorithms were used). The special feature is that the precision needed in this iteration is such 
that any biases are at a level substantially below the statistical noise from the Z2 method. We are able to 
monitor several quantities of interest (eg. TraceA'J -1 on a time slice and the pion propagator to large t) 
continuously during the iterative inversion process. The convergence of these quantities of interest during 
the iterative process is not monotonic, but we are able to establish a value of the residual that guarantees 
sufficiently small systematic errors from lack of convergence. In practice we need approximately one half 
of the number of iterations used in a conventional inversion. This has also been discussed by the SESAM 
collaboration 
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